Can J/tp suppression and pr broadening signal the deconfinement transition at RHIC? 
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We have analyzed the latest NA50 data on J/tp suppression 
in Pb+Pb collisions at CERN SPS. It is shown that a QCD 
based nuclear absorption model, where J ftp's are absorbed in 
nuclear medium could explain the latest NA50 data on the 
centrality dependence of the J/tp over Drell-Yan ratio. The 
model also explains the NA50 data on J/tp over minimum bias 
ratio and thepr broadening of J ftp's. A QGP based threshold 
model where all the J ftp's are suppressed above a threshold 
density, also explains the data sets with smeared threshold 
density. Even at RHIC energy, centrality dependence of J ftp 
suppression or pr broadening could not distinguish between 
the two models. 



I. INTRODUCTION 

J/tp suppression in heavy ion collisions is recognized 
as an important signal of the confinement-deconfinement 
phase transition. NA50 collaboration, at CERN, is a 
dedicated experiment, measuring J/tp cross sections in 
pA/AA collisions. Recently, they have published the 
results of 2000 run of 158 AGeV Pb+Pb collisions [1]. 
The preliminary analysis of 2000 run was presented in 
Quark Matter 2002 [2] . The suppression obtained is still 
anomalous in the sense that the normal nuclear absorp- 
tion model fails to explain it. Compared to 1998 run 
[3], 2000 data are flatter, suppression being more at low 
and intermediate Et- Preliminary analysis indicated less 
suppression (compared to 1998 run) at large Et-, but in 
the final analysis suppression is compatible with the 1998 
run. 1998 NA50 data gave the first indication of anoma- 
lous J /if) suppression and were analyzed in a variety of 
models, with and without the assumption of QGP [4-9]. 
We have shown that a QCD based nuclear absorption 
model, where J ftp's are absorbed in nuclear medium, 
could explain the data [7,8]. We have also shown that 
the model reproduces the NA50 data on the centrality 
dependence of px broadening of J f ip's [10] . What is more 
intriguing is that the predicted J/tp over Drell-Yan ratio 
or the pt broadening at RHIC energy matches with the 
prediction obtained in the QGP based threshold model 
[10]. It seems that even at RHIC energy, Et (centrality) 
dependence of the J ftp suppression or pt broadening may 
not distinguish a deconfining phase transition. 

NA50 collaboration also published the analysis of the 
nuclear absorption of J/tp in high statistics 450 GeV 
pA collisions [11]. They estimated the J/tp nucleon 



absorption cross section (cr^.tf N ) in the framework of 
Glauber model. High statistics 450 GeV pA data yield 
= 4.4±1.0 mb [11]. They also estimate a common 
from latest pA and NA38 200 GeV/c S+U data 
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T abs' N =4.4 ± 0.5 mb. The extracted absorption 
cross section is much smaller than the earlier value of 6.4 
± 0.8 mb extracted from fit to earlier NA50 data [13] 
or 7.1 ± 3.0 mb obtained from a fit to NA38 S+U data 
[12]. Within error, the 200 AGeV S+U cross sections are 
compatible with 450 AGeV pA cross sections. 

In an earlier publication we have analysed the prelimi- 
nary NA50 data of 2000 run [14] . It was shown that QCD 
based nuclear absorption model, with parameters fixed 
from the high statistics pA data give consistent descrip- 
tion of J/tp suppression in 158 AGeV Pb+Pb collisions. 
The preliminary data were also analysed in the QGP 
based threshold model [15]. In the threshold model [4], 
in addition to 'conventional' (Glauber) nuclear absorp- 
tion, an anomalous suppression is used such that all the 
J/tp's are totally suppressed above a critical (threshold) 
density n c . 1998 version of NA50 data [3] were well ex- 
plained in the threshold model, with n c — 3.7-3.75 jmT 2 
and J/f/'-nucleon absorption cross section a^lf N —6A mb 
[4]. J/i/>-nucleon absorption cross scetion 6.4 mb is large 
compared to the recently extracted value of from the high 
statistics pA data, <j'Hf N —AA mb [11]. It was shown that 

with a^lf N =AA mb, the threshold model fails to explain 
the preliminary NA50 data [2] on the centrality depen- 
dence of J/tp over Drell-Yan ratio, unless the threshold 
density is largely smeared [15]. 

In the present paper, we have analysed the latest NA50 
data [1] on J/tp suppression in Pb+Pb collisions . The 
QCD based nuclear absorption model, with parameters 
fixed from the high statistics pA data, still give consis- 
tent description the latest data. QGP based threshold 
model also explains the data if the threshold density is 
smeared. In addition, we have analysed the NA50 data 
on the centrality dependence of J/tp over minimum bias 
ratio [16] and on the centrality dependence of pt broad- 
ening of J ftp's [17]. Both the models, QCD based nuclear 
absorption model and the QGP based threshold model, 
explain these data. The NA50 data of 158 AGeV Pb+Pb 
collisions could not discriminate between the two mod- 
els. It is also shown that even at RHIC energy, centrality 
dependence of J/tp suppression or pt broadening could 
not distinguish between the two models. 

The paper is organised as follows: in section II we 
have briefly described the QCD based nuclear absorp- 
tion model and the threshold model. In section III, NA50 
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data on the centrality dependence of J / tp over Drell-Yan 
ratio and the J/tp over minimum bias ratio are analyzed. 
In section IV, it is shown that both the models could ex- 
plain the NA50 px broadening data. Predicted centrality 
dependence of J/tp over Drell-Yan ratio and pr broad- 
ening at RHIC energy are given in section V. Lastly the 
summary and conclusions are given in section VI. 

II. MODELS FOR J/tp SUPPRESSION 

A. QCD based nuclear absorption model 

In the QCD based nuclear absorption model [7,9], J/tp 
production is assumed to be a two step process, (a) for- 
mation of a cc pair, which is accurately calculable in QCD 
and (b) formation of a J/tp meson from the cc pair, which 
is conveniently parameterized. The J/tp cross section in 
AB collisions, at center of mass energy y/s is written as, 

a J ^(s)=K^ J dq 2 (^p) / dx F <p a/A (x a ,Q 2 ) (1) 

<Pb/B{xb,Q 2 ) Xa * b x F C c^j/^(q 2 ), 

where ^ a b runs over all parton flavors, and Q 2 = q 2 + 
Am 2 . The K factor takes into account the higher order 
corrections. The incoming parton momentum fractions 
are fixed by kinematics and 

W X F + ^Q 2 /* + 

x F )/2 and x b = (y/ 'x 2 F + 4Q 2 /s - x F )/2. a a b^cc 
are the subprocess cross section and are given in [18]. 
Fcc~>j/4)(q 2 ) is the transition probability that a cc pair 
with relative momentum square q 2 evolve into a physical 
J/tp meson. It is parameterized as, 

Fe^j/^q 2 ) = Nj /ip d(q 2 )9(4m /2 - Am 2 c - q 2 ) (2) 

q 2 

(1 ) ap . 

4m' 2 — 4m 2 

In a nucleon- nucleus/nucleus- nucleus collision, the pro- 
duced cc pairs interact with nuclear medium before they 
exit. It is argued [9] that the interaction of a cc pair with 
nuclear environment increases the square of the relative 
momentum between the cc pair. As a result, some of 
the cc pairs can gain enough relative square momentum 
to cross the threshold to become an open charm meson. 
Consequently, the cross section for J / tp production is re- 
duced in comparison with nucleon-nucleon cross section. 
If the J/tp meson travel a distance L, q 2 in the transition 
probability is replaced to q 2 — > q 2 + e 2 L, e 2 being the 
relative square momentum gain per unit length. In [7], 
parameters of the model (a Fl KNj/^ and e 2 ) were fixed 
from experimental data on total J/tp cross section in 
pA/AA collisions, KNj/^0A58, e 2 = 0.225GW 2 / fm 
and a F = 1.0 [7]. 
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FIG. 1. The latest NA50 experimental ratio of the total 
J/tp cross sections and Drell-Yan cross sections in pp and pA 
collisions [11]. The solid lines are the fit obtained in the QCD 
based nuclear absorption model. 

As mentioned in the beginning, recently NA50 collab- 
oration measured J/tp cross section in pA collisions [11]. 
They have measured B tlfJi a(J/tp)/a(DY). In Fig.l, ex- 
perimental data are shown as a function of the nuclear 
length L. The Drell-Yan cross sections donot have any 
A or alternately any L-dependence. The observed L de- 
pendence is then due to J/tp's only. We fit the data 
with two parameters, N norm — KNj ps i/a^Jf(nb) and 
square momentum gain factor e 2 (a F being kept fixed 
at 1). In Fig.l, the fit obtained to the data are shown. 
The two sets of data at 200 GeV/c and 450 GeV/c could 
be fitted with a common square momentum gain fac- 
tor, e 2 = 0.187 GeV 2 /fm, a value 20% lower than the 
value obtained earlier [7]. Lowering of e 2 indicate less 
absorption of J/tp's in nuclear medium, in agreement 
with the Glauber model calculations. While the square 
momentum gain factor do not show energy dependence, 
the evident energy dependence of the cross section ratios 
shows up in the other parameter of the model N norm . 
We obtain N norm = 10.18 at 200 GeV/c and N norm 
= 4.43 at 450 GeV/c. The energy dependence of J/tp 
cross section being taken care of in the model (Eq.l), 
the energy dependence of N norm is due to the Drell-Yan 
cross sections only. In the mass range, 2.9 > M > 4.5 
GeV, the Craigie parameterization [19], of the Drell-Yan 
cross section, a{DY) cx e - 14 - 9M /Vs i gives for the ratio 
a(DY)4 50 GeV /c(£>Y~)2ooGeV = 2.1 -3.1, consistent with 
the presently obtained ratio of 2.29. 

With the parameters of the model fixed from high 
statistic pA data, J/tp production cross section in Pb+Pb 
collisions are obtained following the standard procedure. 
At an impact parameter b, J/tp production cross section, 
as a function of E F is written as [7], 
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J d 2 S T A ( S )T B (b- S )S(L(b,s))P(b,E T ) 

(3) 



where T a .b is the nuclear thickness function, T(b) = 
J_ oo dzp(b, z). P(b,Er) is the Et — b correlation func- 
tion. We have used the Gaussian form for the Et — b 
correlation, 

P(b, E T ) oc exp(-{E T - qN p {b)) 2 /2q 2 aN p (b)) (4) 

where N p (b) is the number of participant nucleons at im- 
pact parameter b. a and q are parameters related to 
dispersion and average transverse energy. For Pb+Pb 
collisions the parameters are, a=1.27 and q— 0.274 GeV 
[4] . S(L) is the suppression factor due to passage through 
a length L in nuclear environment. At an impact param- 
eter b and at point s, the transverse density is calculated 
as, 

n(b, s) = T A (s)[l - e -<^«T B (b- s )] + [A^B] (5) 

and the length L{b,s) that the J/tp meson will traverse 
is obtained as, 



L(b,s) = n(b, s)/2p 



(6) 



At a fixed impact parameter Et fluctuates. Et fluc- 
tuations at a fixed impact parameter plays an important 
role in J/tp suppression. The second drop in the J/tp 
over Drell-Yan ratio, beyond 100 GeV is due to Et fluc- 
tuations [4]. Following Blaizot et al [4], we take into 
account Et fluctuations at a fixed impact parameter b, 
by the replacement: 

L(b, s) -» L(b, s)Et/ <E t > (b). (7) 



B. QGP based threshold model 

We have also analyzed the data in the threshold model 
[4]. The details of the model could be found in [4]. 
Briefly, in the threshold model, in addition to Glauber 
type 'nuclear' absorption, an anomalous suppression fac- 
tor Sanom is used. The J/tp cross section at an impact 
parameter b as a function of Et is then written as, 



d 3 a J /^ 
dE T d 2 b 



J/i> 
'nn 



d 2 sT e J f {s)T e J f {b-s) 
(b,s)P(b,E T ) 



(8) 



where T eff (b) is the effective thickness, T eff (b) 



J x oo dzp(b 7 z)exp(~a a [ s J™ dzlp(b, zl), uj hs is the 
J/V>-nucleon absorption cross-section. In [4], Blaizot ct 



al used a 



J/ipN_ 
abs 



=6.4 mb. In the present analysis, we have 



used o'a^ Ar =4.4 mb, as extracted from the recent pA 



data [11]. In Eq.8, 5' anom (b,s) is the anomalous sup- 
pression factor. Blaizot et al [4] considered two types of 
form for S an om- Assuming that all the J ftp's get sup- 
pressed above a threshold density (n c ), the anomalous 
suppression factor was written as, 



>(b,s) = 0(n(b,s) - n c ) 



(9) 



where n is the transverse density (Eq.5). In ref. [4] it was 
seen that if the theta function is smeared at the expense 
of another parameter, such that suppression is gradual 
rather than abrupt, the quality of fit to data improves 
considerably. This was implemented by writing, 
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( (b, s) = 0.5[1 - tanh(X(n(b, s) - rc c ))] (10) 



In both the form, effect of Et fluctuations at a fixed 
impact parameter was taken into account by rescaling the 
density as, n — > uEt / < Et > (b). The parameters n c 
and A are then obtained by fitting the latest NA50 data 
on centrality dependence of J/ip over Drell-Yan ratio. 

The Drell-Yan pairs do not suffer final state interac- 
tions and the cross section at an impact parameter b as 
a function of Et could be written as, 



d 3 a 



DY 



dE T d 2 b 



- *nn J d 2 sT A (s)T B (b - s)P(b, E T ) (11) 



III. RESULTS 

A. Et dependence of J/tp over Drell-Yan ratio 

In Fig. 2, centrality dependence of J/tp over Drell-Yan 
ratio, as obtained by the NA50 collaboration in their fi- 
nal analysis of the 2000 Pb+Pb run, is shown. Just to 
show, how the Pb+Pb data are changed with time, we 
have also shown the results of 1996-1998 run and the pre- 
liminary analysis of the 2000 run. Final analysis of 2000 
run differ considerably from the first version of the data, 
presumably due to different analysis method. In Fig. 2, 
model predictions for J/tp suppression are also depicted. 

We have used B^^a'^/a^Jj — 38, 17% lower than the 
value obtain from extrapolating 200 AGeV pA/SU data 
(Fig.l) to pp collisions. 200 AGeV pA/SU data is limited 
to L > 3/m and extrapolation to pp data may not be 
very accurate. 

Just to show that the latest data are also anoma- 
lous, we have shown the Glauber model calculation with 
a abs' N '—4.4 mb (the dash-dot-dot line). Only for very 
peripheral collisions, the Glauber model of nuclear ab- 
sorption fits the data. For more central collisions, it pro- 
duces much less suppression than the data exhibit. In 
Fig. 2, the solid line is the calculated ratio in the QCD 
based nuclear absorption mdoel. It agrees well with the 
experiment. The parameters of the model were obtained 
from fitting pA data. In pA collisions we donot expect a 
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deconfining phase transition. Ability of the model to re- 
produce Pb+Pb data, with the same parameters, clearly 
indicate that nuclear absorption alone, treated in an un- 
conventional manner, is capable of explaining the data. 



A 1998 run (final analysis) 

■ 2000 run (preliminary analyisis) 

• 2000 run (final analysis) 

Glauber model 

QCD based nuclear absorption 

Threshold model 

■ ■ Threshold model (smeared) 
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FIG. 2. The experimental centrality dependence of J/tp 
over Drell-Yan ratio in 158 AGeV Pb+Pb collisions. The 
black circles are the final analysis of 2000 run. The pre- 
liminary analysis of the 2000 run is shown as black squares 
and the 1st version of the data obtained in 1996-1998 run 
is shown as black triangles. Different lines in the figure, 
are: (i) Glauber model of nuclear absorption with ^/ iV =4.4 
mb (dash-dot-dot line), (ii) model prediction in the QCD 
based nuclear absorption model (solid line), (iii) one param- 
eter threshold model with threshold density n c =3.78 fmT 2 
(dash-dot line) and (iv) two parameter threshold model cal- 
culation with n c =3.98 fm~ 2 and A=.82 fm 2 (dotted line). 

The QGP based threshold model [4] with only one pa- 
rameter, the threshold density, on the other hand fails 
to give proper description of the data. In Fig. 2 the 
dash-dot-dot line is the best fit obtained to the data 
with threshold density n c =3.78 fm 2 . In the intermediate 
range of Et, agreement with data is not good. Much bet- 
ter fit to data is obtained, when the threshold density is 
smeared. The dash-dot line is the best fit obtained to the 
data with ra c =3.98 /to" 2 and A=0.82 fm 2 . The model 
then reproduces the data through out the Et range. 
Small value of A required for good fit to data indicate 
that considerable smearing of the threshold density is 
required for proper description of the data. The anoma- 
lous suppression is not abrupt but increases gradually 
with density. The threshold density n c we obtain from 
fitting is also larger than the value of 3.7-3.75 obtained 
by Blaizot et al [4]. This is presumably due to smaller 
value of the J/V>-nucleon absorption cross-section. 



B. Et dependence of minimum bias cross section 

In 1996 and 1998 runs, NA50 collaboration obtained 
Et dependence of the J/tp over minimum bias cross sec- 
tions [16]. The minimum bias cross sections are easy 
to calculate. It is essentially the inelastic cross-section. 
In the Glauber model, at impact parameter b and at 
transverse energy Et, the minimum bias cross section is 
written as, 

d 3 MB 

« (1 - exp(-a NN T AB (b))P(b, E T ) (12) 
where T AB (b) = J d 2 sT A {s)T B (s - b). 
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FIG. 3. The experimental centrality dependence of J/tp 
over minimum bias ratio in 158 AGeV Pb+Pb collisions . 
The solid line is the ratio obtained in the QCD based nu- 
clear absorption model. The dash-dotted line is the threshold 
model calculation with n c =3.98 fmT 2 and A=.82 fm 2 . 

In Fig. 3, the NA50 data [16] are shown. Data are 
given in arbitrary units. The QCD based nuclear ab- 
sorption model with parameters fixed from pA/AA total 
J/tp cross section (the solid line) well describe the data. 
Only at large Et the predictions differ from the experi- 
ment. Data at very large Et may not be reliable. Final 
analysis of 2000 run show less suppression of J/tp than 
found in the analysis of 1996-1998 run (see Fig. 2). The 
J/tp over minimum bias ratio is also well explained in the 
two parameter threshold model (the dashed line). Like 
the J/tp over Drell-Yan ratio, the J/tp over minimum bias 
ratio is also incapable of distinguishing between the two 
models. 



IV. P T BROADENING OF J/tp 



It is well known that in pA and AA collisions, the sec- 
ondary hadrons generally shows &pt broadening [20,21]. 
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Kharzeev ct al [21] suggested px broadening as a probe 
for the deconfining transition. They argued that in a de- 
confining medium, px broadening will visibly decrease at 
large Et, in contrast to a gradually increasing pt broad- 
ening in a nuclear medium. Recently we have shown 
that decreasing pt broadening can not be considered as 
a signal of deconfining phase transition, as such a trend is 
also obtained in the QCD based nuclear absorption model 
[10]. It was also shown that Et fluctuations, at a fixed 
impact parameter plays an important role in explaining 
the experimental data on p T broadening of J/ip's. The 
QCD based nuclear absorption model [7] could explain 
the NA50 pr broadening data if the effect of Et fluctua- 
tions are properly accounted for. In [10], we have also ob- 
tained pt broadening in the threshold model, with a sin- 
gle parameter, the threshold density, n c =3.7 fm~ 2 and 
a abt N '=6-4 mb. One parameter threshold model could 
not fit the data. As it was shown in last two sections, 
with the latest NA50 data [2,11], parameters of both the 
models are changed and it interesting to see the conse- 
quence on pt broadening. 

The natural basis for the pt broadening is the initial 
state parton scatterings. For J/^'s, gluon fusion being 
the dominant mechanism for cc production, initial state 
scattering of the projectile/target gluons with the tar- 
get/projectile nucleons causes the intrinsic momentum 
broadening of the gluons, which is reflected in the pt 
distribution of the resulting J/^'s. Parameterizing the 
intrinsic transverse momentum of a gluon, inside a nu- 
cleoli as, 



f(q T ) ~ exp(-q T / < q\ >) 



(13) 



momentum distribution of the resulting J/ip in NN col- 
lision is obtained by convoluting two such distributions, 



exp(-p T / < P T >nn) 



(14) 



where < p% > J jl%=< q T > + < >■ In nucleus- 



q T 



nucleus collisions at impact parameter b, if before fu- 
sion, a gluon undergo random walk and suffer N num- 
ber of subcollisions, its square momentum will increase 
to q T — > q T + NSq, S being the average broadening in 
each subcollisions. Square momentum of J/ip then easily 
obtained as, 

< Pt >ab (b) =< p 2 t >Ht +S N AB (b) (15) 

where Nab (b) is the number of subcollisions suffered by 
the projectile and target gluons with the target and pro- 
jectile nucleons respectively. 

Average number of collisions Nab (b) can be obtained 
in a Glauber model [21]. At impact parameter b, the po- 
sitions (s, z) and (b — s, z') specifies the formation point 
of cc in the two nuclei, with s in the transverse plane and 
z, z' along the beam axis. The number of collisions, prior 
to cc pair formation, can be written as, 



N(b, s, z, z') = a gN / dz A pA{s,z A ) (16) 

J —CO 

+cr gN / dz B p B (b - s, z') 

J — oo 

where a g N is the gluon-nucleon cross section. Above ex- 
pression should be averaged over all positions of cc for- 
mation with a weight given by the product of nuclear 
densities and survival probabilities S, 

//•OO /• oo 

d 2 s dz PA {s,z) dz' PB (b-s,z')x (17) 
J — oo J —oo 

S(b,s,z,z')N(b,s,z,z')/ / d 2 s / dzp A (s,z)x 

J J —oo 

/oo 
dz'psib — s, z')S(b, s, z, z') 
-oo 

Finally, corresponding quantity at fixed transverse en- 
ergy Et is obtained as, 

Nab(E t )= J d 2 bP(b,E T )(TABNAB(b)/ (18) 

j d 2 bP(b,E T )<7AB 

where <jab is the inelastic cross section for AB collisions. 
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FIG. 4. (a) NA38 experimental data on the centrality de- 
pendence of J/ip over Drell-Yan ratio, in 200 GeV/c S+U 
collisions. The solid line is the ratio obtained in the QCD 
based nuclear absorption model, (b) The centrality depen- 
dence of pr broadening in S+U collisions. The solid line is a 
fit to the data in the QCD based nuclear absorption model. 
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Fluctuations of Et at a fixed impact parameter will 
also affect the average number of collisions Nab{Et)- 
We have taken into account the Et fluctuations by the 
replacement, 

N AB (b) -> Et/ <E t > (b)N AB (b). (19) 
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FIG. 5. The energy dependence of < p% >nn along with 
the fit to it with Eq.20. 

Pt broadening of J ftp's in A A collisions depends on 
two parameters, (i) < p T >nn, the mean squared trans- 
verse momentum in NN collisions, a measurable quantity 
and (ii) the product of the gluon-nucleon cross section 
and the average parton momentum broadening per colli- 
sion, dgNSo- Since gluons are not free, the second quan- 
tity is essentially non measurable. We obtain a g N5o from 
a fit to the NA38 pt broadening data [12] in S+U colli- 
sions at 200 GeV/c. < p T >nn at corresponding energy 
is known from NA3 experiment, < p\ >nn— 1.23 ±0.05 
[22]. The Et — b correlation parameters, a and q for 
S+U collisions are, a=3.2 and g=0.74 GeV [23]. To 
show that the present 'unconventional' nuclear absorp- 
tion model also reproduces the centrality dependence of 
J/tp over Drell-Yan ratio in S+U collisions, in Fig. 4a, we 
have compared our results with the experimental data. 
We have neglected the effect of Et fluctuations in S+U 
collisions. The agreement between data and theory is 
good. In Fig. 4b, NA38 experimental data on the Et de- 
pendence of pt broadening are shown. The solid line is 
a fit to the data, obtained with < p\ >nn— 1-23 (fixed) 
and (jgN&o = 0.442 + 0.056. Value of dgNSo agrees closely 
with the value obtained by Kharzeev et al [21] in the 
conventional nuclear absorption model and also with the 
value obtained in the comover model [24] . 

< Pt >nn increases weakly with energy. To obtain 
< Pt >nn for Pb+Pb collisions at 158 GeV/c, we have 
fitted the existing experimental data [22,25-27] with log- 
arithmic energy dependence, 

< P T >nn= a + b\n^/s (20) 



In Fig. 5, experimental data along with the fitted curve 
obtained with a = —0.38 and b = 0.53. is shown. From 
the above parameterization, we obtain, < p T >nn =1.15 
GeV 2 , for Pb+Pb collisions at CERN SPS. As we intend 
to predict pt broadening at RHIC energy, < p T >nn at 
RHIC energy (-^3=200 GeV) is also obtained from the 
above parameterization. At RHIC energy, < p\ >nn 
=2.45 GeV 2 . However, we must warn our reader to treat 
the above number with caution. The experimental data 
being limited to 60 GeV only, extrapolation to RHIC 
energy is unreliable. 
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FIG. 6. The experimental centrality dependence of pr 
broadening in Pb+Pb collisions. The solid and dashed lines 
are the QCD based nuclear absorption model predictions, 
with and without the effect of Et fluctuations on the average 
number of gluon-nucleon collisions JVab(&)- The dash-dot-dot 
and dash-dot lines are the same in the QGP based one param- 
eter threshold model. The prediction in the two parameter 
threshold model is shown as dotted line. 

In Fig. 6, we have shown the result of pt broadening in 
the model. The solid and dashed lines are the < p\ > 
in our QCD based nuclear absorption model, with and 
without the effect of Et fluctuations on Nab (b) ■ When 
the Et fluctuations are not accounted for (the dashed 
line), the model could not explain the experiment be- 
yond Et=100 GeV. Experimentally, < p T > continues 
to increase with Et beyond 100 GeV, but the model pre- 
dicts a decreasing trend. Beyond 100 GeV (the knee 
of the Et distribution), J ftp's are strongly suppressed. 
Strong suppression causes the pt broadening to decrease 
beyond 100 GeV. The decreasing trend is changed into an 
increasing trend if the effect of Et fluctuations on Nab is 
taken into account (the solid line) . Et fluctuations effec- 
tively increases the average number of collisions NabQ)) 
and counter balance the strong suppression effect beyond 
the knee of the Et distribution. Considering that all the 
parameters of the model were fixed, model describe the 
data very well. 

In Fig. 6, the one parameter threshold model prediction 
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with threshold density n c =3.78 jmT 2 is shown. When 
the effect of Et fluctuations on Nab is neglected (the 
dash-dot-dot line) the model fails to explain the data. 
At low Et, it predict < p T > in accordance to the QCD 
based nuclear absorption model, but beyond Et=60 
GeV, it predict less pt broadening. Also, the increas- 
ing tendency beyond 100 GeV is not reproduced. Even 
when the effect of Et fluctuations are taken into account 
(the dash-dot line), the model fails to give proper de- 
scription to the data. This is expected as the one param- 
eter threshold model donot give very good description to 
the centrality dependence of J/ip over Drell-Yan ratio. 
Hucfner et al [28] also analyzed the NA50 pt broadening 
data in the threshold model and essentially obtain a simi- 
lar result that the model could not fit the NA50 pt broad- 
ening data. As shown earlier, a two parameter threshold 
model, with smeared threshold density, reproduced the 
centrality dependence of J/ip over Drell-Yan ratio as well 
as J/ip over minimum bias ratio. In Fig. 6, the prediction 
obtained in the two parameter threshold model, with the 
effect of E T fluctuations included, is shown as the dotted 
line. Considering that it is also a parameter free calcula- 
tion (< p T >nn and (JgN^a fixed), the model describes 
the data well. Here, we may mention that Armesto et 
al [24], in the comover model also explained the NA50 
Pt broadening data. However, as they did not account 
for the Et fluctuations on NAs{b) the increasing trend 
beyond 100 GeV could not be reproduced. 



V. PREDICTION FOR RHIC AU+AU 
COLLISIONS 

Recently PHENIX collaboration published the central- 
ity dependence of charm production in Au+Au collisions 
at RHIC energy [30] . Centrality dependence of charm 
quark production is consistent with binary collisions scal- 
ing. PHENIX collaboration also published the yield of 
J/ip in a few centrality ranges of Au+Au collisions at 
RHIC energy [31]. Data have very large error bars. Data 
do not show any indication of large enhancement as spec- 
ulated in some models [32]. We have shown that the 
PHENIX data on the centrality dependence of J/ip pro- 
duction are well described in the QCD based nuclear ab- 
sorption model [33]. In this section we compare the J/ip 
production at RHIC energy in the QCD based nuclear 
absorption model with the production in the QGP based 
threshold model. Model parameters are kept fixed at 
the values required for Pb+Pb collisions at SPS energy. 
At RHIC energy the so called hard scattering, propor- 
tional to number of binary collisions appear. However, 
PHENIX data on J/ip production [31] do not require any 
hard scattering component [33] and presently we neglect 
them. 

In Fig. 7, we have compared the suppression factor for 
J/ip production at RHIC energy in the QCD based nu- 
clear absorption model and in the QGP based threshold 
model. In a centrality range of collisions, J/tp suppres- 
sion factor is defined as, 
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FIG. 7. Suppression factor for J/ip production as a func- 
tion of average number of NN collisions in Au+Au collisions 
at RHIC energy. The solid and dashed lines are the predic- 
tion from QCD based nuclear absorption model and the QGP 
based threshold model respectively. In the figure, the filled 
circles represent the PHENIX data on centrality dependence 
of total charm production (scaled by a factor of 1500) per NN 
collisions. 



S = 



J/4> 

AA 



< N AA > a J J% 



(21) 



where < N >aa is the average number of NN collisions 
in the centrality range of collisions. As seen in Fig. 7, 
both the QCD based nuclear absorption model and the 
QGP based threshold model predict nearly same sup- 
pression factor. In Fig. 7 we have also shown the recent 
PHENIX measurement of total charm quark multiplic- 
ity per NN collisions (N cS /Tab) [30]. Data were scaled 
by a factor of 1500. Except for the very peripheral colli- 
sions (60-92% centrality), centrality dependence of charm 
production agree well with the dependence predicted in 
both the models. Possibly centrality dependence of J/ip 
production will not distinguish between the models. 

As told earlier, pt broadening of J/ip is considered as a 
probe of deconfinement transition. We have seen that at 
SPS pt broadening do not distinguish between the QCD 
based nuclear absorption model and QGP based thresh- 
old model. In Fig. 8, model predictions for p T broadening 
of J/ij/s at RHIC energy are shown. Both the models 
predict very similar pt broadening. As it is in SPS en- 
ergy, at RHIC energy also pt broadening of J / ip do not 
distinguish between the QCD based nuclear absorption 
and QGP bases threshold models. 
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FIG. 8. pr broadening of J/ip's as a function of average 
number of NN collisions at RHIC energy. The solid and 
dashed lines are the prediction from QCD based nuclear ab- 
sorption model and the QGP based threshold model respec- 
tively. 



VI. SUMMARY AND CONCLUSIONS 

To summarize, we have analyzed the latest (2000 run) 
NA50 data on J/tp suppression in 158 AGeV Pb+Pb col- 
lisions. QCD based nuclear absorption model, with pa- 
rameters fixed from the NA50 high statistics pA data, 
well explain the latest NA50 data. The model also ex- 
plain the centrality dependence of J/tp over minimum 
bias ratio as well as the centrality dependence of px 
broadening of J ftp's. The same data sets were analysed in 
the QGP based threshold model. Threshold model with 
smeared threshold density also explain those data sets. 
At SPS energy J ftp's can not distinguish between the two 
differeny models. We have also shown that even at RHIC 
energy, both the models predict nearly similar centrality 
dependence of J/tp suppression and its pr broadening. 
Possibly, for the deconfinement phase transition, there is 
no 'smoking gun'. A variety of data, analysed in a con- 
sistent manner will be able to shed light on the possible 
deconfinement phase transition. 
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